Enhanced high-energy neutrino emission from choked gamma-ray bursts due to meson 
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It has been suggested that a potentiaUy large haction of supernovae could be accompanied by 
relativistic outflows that stall below the stellar surface. In this letter we point out that internal 
shocks that are believed to accelerate protons to very high energies in these flows will also accelerate 
secondary mesons and muons. As a result the neutrino spectrum from meson and muon decay is 
expected to be much harder compared to previous estimates, extending as a single power law up to 
~10^ TeV. This greatly improves the detection prospects. 
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Based on the observational connection (see Ref. [l| 
for a review) between gamma-ray bursts (GRBs) and su- 
pernovae (SNe), it has been hypothesized that a sizable 
fraction of all SNe is accompanied by a relativistic out- 
flow similar to those that are believed to lie at the base of 
observed GRBs ^] . Initially the flow accelerates by radi- 
ation pressure and ploughs through the pre-burst stellar 
material. The flow may however be 'choked' below the 
stellar surface when the central engine is not active for a 
sufficiently long time. As electromagnetic radiation that 
is dissipated by these flows will be absorbed by the stel- 
lar material, neutrinos are likely the only particles that 
could indicate the existence of this phenomenon. 

High-energy (> 1 TeV) neutrinos arise predominantly 
in the decay of charged mesons and muons that are cre- 
ated in the interactions of shock-accelerated protons with 
target protons or photons. In particular, neutrino emis- 
sion due to proton acceleration in internal shocks in the 
relativistic flow has been studied in detail [H, 13, 0, S 0] ■ 
It has been argued that the fluence of high-energy neu- 
trinos in this scenario is suppressed because the mesons 
and muons lose a substantial amount of energy before 
decay due to synchrotron emission and hadronic interac- 
tions 0, 0| . This energy loss strongly limits the neutrino 
detection prospects. However, as we point out here, the 
mesons and muons are also deflected and scattered by 
the strong magnetic field on very short timescales. These 
particles have ample time to cross the shock several times 
before decay and hence they are subject to shock accel- 
eration. In this work we investigate the observational 
consequences of the acceleration of neutrino parent par- 
ticles. For concreteness we restrict ourselves to internal 
shocks in choked GRB outflows. We expect however that 
the acceleration mechanism is quite general in GRBs, so 
that our results may affect estimates of neutrino fluxes 
in other scenarios as well. 

The model — We consider a relativistic outflow with 
total energy E ~ 10^^i?52, Lorentz factor Tj — lOFj,! 



and opening angle 9 = O.10-i (we use the notation 
Q = l^fQx throughout this letter). Here parame- 
terizes the combined effect of collimation and relativis- 
tic beaming. The isotropic-equivalent burst energy is 
= ^E/e"^ = 2 X 10^'' erg x Ozl £^52- Following Refs. 
[j, d, we assume that internal shocks occur in the 
flow at a radius Tint — 2cF|^t = 6 x 10^^ cm x r| <5t-i, 
where 5t = 0.1(5j^_iS is the variability timescale of the 
central engine. The comoving proton density at the 
internal shock radius is n'^ = i?iso/47rr?jjF|rnpC^i = 
10^^ cm~3 X Yjl ezl E52 t^^ Stzj, where t = Wti s is the 
burst duration (quantities in the comoving frame are de- 
noted with a prime) and rup is the proton mass. The 
large proton density gives rise to a very high Thom- 
son optical depth r > 10^. This implies that syn- 
chrotron photons that are emitted by shock-accelerated 
electrons will thcrmalize. The number density of photons 
at the internal shock radius is n'^ — Q-HiteU' jhcfl^ — 

2 X 10^3 cm-3 X FT;'/^ ^^3/2 (,^^)3/4^^3/4 ^j^^^^ 

U' = i?iso/47rr?jjF|ci is the comoving energy density and 
Ce = O.lee,-! denotes the fraction of the total energy in 
the flow that is transfered to the thermal photon distribu- 
tion. The magnetic field strength at the internal shock 
radius is B' ^ 2 x 10* G x T-lez\{tBE)y2^t~^''^ 5tz\, 
where es — O.les.-i denotes the ratio of electromagnetic 
energy to the total energy in the flow. 

We adopt the SN rate parameterization presented in 
Ref. The SN rate within proper distance dp can be 
approximated with iVgN 4 x 10^ year~^ ^ 21 where 
dp = 100 dp, 2 Mpc. This estimate is within ~30% for 
dp < 500 Mpc. Due to collimation of the relativistic out- 
flow, the rate of observable choked GRBs associated with 
these SNe is TVcgrb = O^Nm/'^ = 2 yr'^ x Csn^' id^,2, 
where ^sn < 1 is the fraction of SNe that is endowed with 
the type of outflows considered in this work. 

Proton acceleration — We assume that internal 
shocks accelerate a fraction of the protons in the flow 



2 



to high energies (see also below). We express the en- 
ergy spectrum of accelerated protons as dNp/dE^ = 
^pE,,o{p-l)TJ^E'p-P{mpC^)P-^, where = O.Ol^p^^s 
denotes the fraction of shock-accelerated protons to all 
nucleons in the flow, and p is the shock-acceleration 
power-law index. Theoretical studies Q indicate that 
p ~ 2.3, but recent observations suggest that p could 
vary To keep the discussion general we consider 

here the range p = 2.0 . . . 2.6. The proton acceleration 
timescale is equal to ip — ^'pIl^B'T'^ = 6 x 10^^^ s x 

E'j,,,T'l^e_i{tBE)-^'^5t-it\'^, where E'^^ = lE'^^^TeW 
denotes the proton energy and = lOF'^, i is the Lorentz 
factor of the shock Q. We take — Tj since the vari- 
ation in Lorentz factors between two subsequent shells 
of material AFj Tj. The maximum proton energy 
may be limited both by energy losses and by the fi- 
nite shock lifetime. The lifetime is roughly equal to the 
dynamical timescale t'^y^^ = Hnt/cFj = 2s x Tj^i5t-i. 
The dominant proton energy-loss mechanisms are syn- 
chrotron radiation, photopion production and proton- 
proton [pp) collisions. The synchrotron energy-loss 
timescale is sync — ['oT^rn^c^) / (aTmlB''^ E'^) = 0.1s x 
^'p}i^%i^^-i^^BE)'^l5t^_iti, where we assume that the 
protons are relativistic. The energy-loss timescale due 
to pp coUisions is ip,pp = \/{cKppapp{l — £,p)n'p) — 
10-^ s^Tl^e'i^E^^6t^_^ti, where we assume that ^p ^ 1 , 
and we approximate the cross section for pp collisions 
with (Tpp ~ 5 X 10^^^ cm^ and the fractional energy loss 
with Kpp ~ 0.5. At center-of-mass energies well above 
the pion creation threshold we approximate the proton- 
photon (p7) cross section with Cp-y ~ 10~^^ cm^, and the 
fractional energy loss to pion production with Kpjj^ ~ 
0.2. In this regime the bulk of the photons participates 
in photopion production so that we estimate the photo- 
pion energy-loss timescale as tp p^^ ~ 1 / (cap-^Kp-^-j^n'^) = 

9 X 10-6 s X Y'^l^&^ll{t^E)lll^8f-llt'\'^ . Equating the ac- 
celeration timescale to the dynamical timescale and the 
energy-loss timescales we find that, unless extreme val- 
ues of the parameters are invoked, the maximum proton 
energy E'^ ^^^^ is determined by synchrotron energy loss. 
In this case 

K,^.. = 5 X lO^TeV X F2^0V2(,5iJ)-V45iV2^i/4 

For protons with sufficient energy the optical depth 
for photopion production Tp^ ~ 1ap^n!^r\ntlV j = 2 x 

10^ X Y^y'^rJ^'^{teE)ll^Stl\''^tfl'^ is always larger 
than unity (unless Fj is very large). On the other 
hand, the optical depth for pp interactions Tpp ~ 
CTppn^rint/Fj = 3 X 10"^ X V-l6z\E^25t~_\t'{^ is less than 
unity when Vj ^80. In this case protons with energy 
below the photopion threshold may traverse the jet rel- 
atively unhindered and impact directly on the jet head 
or the stellar material, which could have interesting ob- 
servational consequences. Here we assume that protons 



lose all their energy in the outflow and that 20% of this 
energy is transfered to secondary mesons. 

Meson and muon acceleration — The secondary 
mesons and muons are created in a strongly magne- 
tized environment. The magnetic field deflects and scat- 
ters the particles through electromagnetic interactions on 
timescales comparable to the gyration timescale t'^ gyj. ~ 

eUqcB' = 6 X lO-^^s x e^^pF^ ^6'_l(eB£;)^//^(5^_l^y^ 
where x denotes either fj, (muon), tt (pion) or K (kaon), 
and = 1 p TeV is the particle energy. The co- 
moving decay time is given by t'^ = T^f'^j-mxC? = 
2 X 10^^ s (2 X lO""' s, 2 X 10^^ s) x g for muons (pions, 
kaons), where t^^ denotes the proper decay time and mx 
denotes the mass. Since the decay time is much longer 
than the gyration time, the particles have ample time to 
be deflected and scattered by the magnetic held. This al- 
lows them to cross the shock repeatedly, thereby gaining 
energy in a stochastic way. This is essentially the same 
mechanism of shock acceleration that applies to protons. 
A necessary condition for this mechanism to work is that 
the acceleration timescale ^cc = ^'x.^jxl^'s ^ dec- 
Notice that the ratio t'x.acd^'x dec independent of the 
particle energy, so that there is no intrinsic maximum 
energy to the acceleration process. 

For stable particles the energy spectrum due to shock 
acceleration can be approximated with a power law with 
index p = 1 — ln(7'i-et)/lnx, where V-^ct — 0.5 is the return 
probability (i.e. the probability that a particle completes 
a full cycle of two shock crossings), and x — 1-6 — 2.0 is 
the average fractional energy gain per cycle [8^. For un- 
stable particles the energy spectrum at decay (which de- 
termines the energy spectrum of the daughter particles) 
is, in principle, expected to be softer than the energy 
spectrum of accelerated stable particles because fewer 
particles complete many cycles. The effect of particle 
decay can be accounted for by adopting the return prob- 
ability Vrct — Vrct — ^dec, whcrc V-cd IS the return prob- 
ability for stable particles and "Pdec — t'x acc/^L dee 
probability that a particle decays during one cycle. The 
acceleration timescale of charged particles in the internal 
shock environment is <^_aee = ^UqcBT'^ = 6 x 10"^-^ s x 
<,o^liO-i{eBE)-^'/^5t.itY\ Since t^,, » t^,, over a 
wide range of parameters, "Pdec ^ "Pret and we expect 
that the mesons and muons are accelerated to a power 
law with index p ~ 2.3 that extends to the maximum en- 
ergy determined by the shock lifetime or by energy losses. 
We find that, both for mesons and muons, the maximum 
energy e'^^^^ is determined by synchrotron losses. Equat- 
ing the acceleration timescale to the synchrotron energy- 
loss timescale i^.syne = tp.synd^x/'mp)'^, we estimate: 

e;'""'^ = 6 X 10^ TcV X Vl^e"/^ {eBE)-,'/'St'/^tY' ; (2a) 
^^max ^ ^Qs^ev X rl,0'/^{eBE)-,'/^5t'/^tY';{2h) 
= lO^TeV x Vl^e"/^ ieBE)-,'^^6t'/^tY\i2c) 
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Neutrino fluence — We denote the average neutrino 
multiplicity per proton with M.p^(^x) and the average frac- 
tion of the proton energy that is transferred to the neu- 
trino with . Here x labels the intermediate particle 
(muon, pion, kaon). Following Refs. @, 0], we assume 
that a proton transfers 20% of its energy per collision to 
the secondary mesons. Furthermore we take the average 
pion (kaon) multiplicity per proton interaction equal to 
1 (0.10). The branching ratio of pion (kaon) decay to a 
muon and a neutrino is virtually unity (0.63) and the neu- 
trino receives ~0.25 (0.50) of the meson energy. A muon 
transfers ~0.33 of its energy to each of two daughter neu- 
trinos. Hence Mp^(f,) = 2, Mpy^^^) = 1, Mp^(K) = 0.06; 

^pi^in) = 0.05, ^pi/(7r) = 0.05, Cpiy{K) =0.1. 

The differential neutrino fluence in the observer frame 
(all flavours combined; neutrinos and antincutrinos com- 
bined) can be expressed as follows: 



1 1 Mp^(x) dNp 



A-Kd^Tj ^p^(^) dEp 

^u{x)'r^j^l^&ZlE52£.padp,2 



(3) 



where we assume that the redshift z ^ 1, and 

$^(^) = 5.2 X lO^'^TeV^^cm-^ (4) 

From eqs. ([2]) , we find that the maximum neutrino energy 
in the observer frame is: 



^max 



- 2 X 10^ TcV X Tl,e'/^ieBE,,)-'/''St'/^tl/^ ; 
e-(7) = 3 X 10^ TeV x Vle'J^ {esE^.r^/^JX^' ; 
e',7;J) = 6 X lO^TeV x Ty/,\eBE,,)-^^^St'/^t\/' , 

where we take the neutrinos to be isotropic in the comov- 
ing frame. 

Detection prospects — Based on preliminary results 
presented in Ref. flo| we conservatively approximate the 
effective area of IceCube for muon neutrinos with: 



^cff(eiy) 



5.0 X 102 cm2 X (e.,o)'-^ (e?h < < ^br) 



1.5 X lO^cm^ X (e^,o)°-^ (e^ > ebr) 



where Cth = 0.1 TeV denotes the detector threshold en- 
ergy, ebr — 60 TeV is a break energy, and = 1 e„,o TeV. 
As the neutrino flavour ratio at the source is roughly 
I'e : i'^ : i^r — 1 : 2 : 0, the expected neutrino flavour ratio 
at the detector is ~ 1 : 1 : 1 due to neutrino oscillations 



over very large distances llj . Hence we approximate the 
fluence of muon neutrinos at the detector with 

{Tr)i<^i^) + ^i'{K){<^i^)) ■ (5) 

We estimate the number of muon-neutrino interactions 
Ni, in IceCube by multiplying the muon-neutrino fluence 



with the effective area and integrating over to find: 

N,^ ^.,N,TPf9zlE,2^p,2d;'2 , (6) 
where = ($^(^) + + *i.(K))/3, and 

9 X 10^ TeV cm^ / 1 ~ x^Y"^ 1 - xldz^ \ 



60P- 



2.7-p 1.3-p 



Here we use the shorthand notation xth = eth/cbr and 
Xmax = e™'''^/ebr, and assume that e^'^^ > ebr and that 
1.3 < p < 2.7. In the limit that e™^'' > ebr, we find that 
N^^ = 2 X 10^ (8 X 10^, 4 X 10^) TeV cm^ for p = 2.0 
(2.3, 2.6). Combining this with eqs. ([4]) and ([61), we pre- 
dict N^^ = 14 (0.63, 0.038) V'^^^^ezlE^2^p,^2dp^ muon- 
neutrino interactions in IceCube for p = 2.0 (2.3, 2.6). 
Hence, for model parameters similar to those adopted in 
this study, a choked GRB at 100 Mpc could be observed 
by IceCube provided that the shock-acceleration index p 
is not too large. For reference values of the other param- 
eters the detection of one neutrino requires p < 2.3. 

The diffuse flux (per sterad) due to unresolved choked 
GRBs can be estimated with 
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87r 



dz 



dV_ 

dz 



nsn^u^i^v). (7) 



where ^sn denotes the fraction of SNe that is accompa- 
nied by a choked GRB, V is the comoving volume, z is 
the redshift, and nsN in the SN rate per unit volume for 
which we adopt the parameterization given in Ref. Q 
(see also Ref. @). We find that: 



<(e.) = |.^fr^~2i?5,Cp.2&N {-^) 



(8) 



where $5^^^ = 9 x 10^^ $1,^, sr^^. In figure[T]we have plot- 
ted our estimate ([8]) of the diffuse muon-neutrino flux for 
three values of p and reference values of the other pa- 
rameters. Also shown are the 90% confidence level upper 
limits of the AMANDA-II [l?| and IceCube (3 year) [l|| 
experiments, the atmospheric neutrino background with 
the parameterization used in Ref. and the Waxman- 



Bahcall bound [14[. As can be seen in the figure, the ex- 



isting limit from the AMANDA-II experiment is already 
constraining the parameter space of choked GRBs. For 
p = 2.0 (2.3), we find that Er^2S,p.2^SN S 10"^ (!)■ In this 
regime the predicted diffuse flux is above the Waxman- 
Bahcall bound, which applies to optically thin sources. 
When p = 2.6 the expected diffuse flux is below the 
detector sensitivity for reference values of the other pa- 
rameters. IceCube will be able to put more stringent 
constraints on the model parameters. A more detailed 
analysis of the detection prospects is beyond the scope 
of this work. We note however that a visible SN coun- 
terpart, which would provide evidence in favour of our 
model, would strongly reduce the neutrino background. 
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FIG. 1: Diffuse muon-neutrino flux at Earth for three 
values of the power-law index p and reference values of the 
other parameters. Also shown are the atmospheric neutrino 
background, the AMANDA-II (A) and IceCube (I) upper 
limits, and the Waxman-Bahcall bound (WB). 

It may also be feasible to use neutrino detectors to initi- 
ate a SN search. 

Conclusions — In this letter we have found that 
muons, pions, and kaons created by pp and P7 interac- 
tions in the internal shocks of choked GRB outflows are 
also accelerated by these shocks, essentially in the same 
way as protons and electrons are. Although this accel- 
eration mechanism appears to be unavoidable under the 
circumstances believed to be present in these flows, it 
was (to the best of our knowledge) not considered be- 
fore. Due to meson and muon acceleration the resulting 
neutrino spectrum is expected to be dominated by neu- 
trinos from muon decay and to follow a single power law 
with index p ~ 2.3 up to a maximum energy ^10'^ TeV 
(and larger than 10** TeV for the subdominant contribu- 
tion of neutrinos from kaon decay). This is in contrast to 
previous estimates that predict spectral breaks at ener- 
gies < 1 TeV for neutrinos from meson decay and virtu- 
ally no high-energy neutrinos from muon decay 0, H, [1] . 
The relatively hard neutrino spectrum strongly increases 
the detection prospects. In fact, the current AMANDA- 
II limit on the diffuse neutrino background is already 
mildly constraining the model parameters. The upcom- 
ing IceCube neutrino detector will be in a good position 
to test the model and, possibly, to observe neutrino emis- 
sion from single choked GRBs. We have estimated that 
a single choked GRB at 100 Mpc aimed toward Earth 
will result in 14 (0.63, 0.038) muon-neutrino interactions 
in IceCube for p = 2.0 (2.3, 2.6) and reference values of 
the other parameters. The rate of choked GRBs emitting 
neutrinos toward Earth within 100 Mpc may be as large 
as a few per year. 

An important caveat in our results is that the model 
relies on the existence of internal shocks that occur at a 
radius Tint due to variability in the flow. It is however 
not clear whether these shocks can indeed develop while 



the jet is traversing the star and has not yet created a 
low-density funnel (Thomas Janka, private communica- 
tion). In contrast to this, a forward shock and a reverse 
shock seem unavoidable in the interaction of the relativis- 
tic outflow with the stellar environment. Also internal 
shocks are expected to occur behind the forward shock 
as it propagates through the star. Neutrino production 
in these shocks will be studied in a forthcoming publica- 
tion. We have also assumed that a fair fraction of the 
secondary mesons makes its way to the shock after be- 
ing produced. This has to be verifled in a more detailed 
study. We expect that this assumption is best justifled 
at high proton energies, where the proton mean free path 
is not much larger than the meson gyroradius. 

Shock acceleration of mesons and muons may be a 
fairly general phenomenon in GRBs. The necessary con- 
dition that the acceleration timescale is smaller than the 
decay timescale is fulfilled when BTj > 5 (6x10^, 5 x 10^^) 
G for muons (pions, kaons), which is easily achieved in 
GRBs. Hence our results may also affect the estimates 
for neutrino emission from successful GRBs. 
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